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Abstract 

In this paper we briefly present our works on the relic gravitational waves (RGW) 
and the CMB polarization in the accelerating universe. The spectrum of RGW has 
been obtained, showing the influence of the dark energy. Compared with those 
from non-accelerating models, the shape of the spectrum is approximately similar, 
nevertheless, the amplitude of RGW now acquires a suppressing factor of the ratio of 
matter over dark energy oc Q m /Q\ ~ 0.4 over almost the whole range of frequencies. 
The RGW spectrum is then used as the source to calculate the spectra of CMB 
polarization. By a two half Gaussian function as an approximation to the visibility 
function during the photon decoupling, both the "electric" and "magnetic" spectra 
have been analytically derived, which are quite close to the numerical ones. Several 
physical elements that affect the spectra have been examined, such as the decoupling 
process, the inflation, the dark energy, the baryons, etc. 
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1 Introduction 

The existence of gravitational waves is a major prediction of General Relativity that has not 
yet been directly detected. On other hand, inflationary models predict, among other things, a 
stochastic background of relic gravitational waves (RGW) generated during the very early stage 
of expanding universe. [1]~[7] Therefore, the detection of RGW plays a double role in relativity 
and cosmology. For a number of gravitational detections, ongoing or under development, 
the spectrum of RGW represents one of their major scientific goals. However, the current 



expansion of the universe has been found to be an accelerating one, probably driven by dark 
energy. This will have important implications on RGW and its detections. As is known, 
the cosmic background radiation has certain degree of polarization generated via Thompson 
scattering during the decoupling in the early universe. [8] [11] In particular, if the tensorial 
perturbations (RGW) are present at the photon decoupling in the universe, then magnetic 
type of polarization will be produced. [12]~ [25] This would be a characteristic feature of RGW 
on very large scales, since the density perturbations will not generate this magnetic type of 
polarization. Besides the generation of linear polarization, the rotation of linearly polarized 
EM propagation by RGW has also been first studied in Refs.[26, 27], and WMAP polarization 
data has already been used to constrain the effect in cosmological distance to 0.1 rad, which 
is important for fundamental physics. For both theoretical and observational studies, it is 
necessary to examine the effects of the dark energy on RGW and on CMB anisotropies and 
polarization. In this talk I shall present our calculational results on these issues. 

First I will present briefly our result of the spectrum of RGW, both analytical and numerical, 
in the accelerating Universe f2 A + Vt m = 1. As a double check, we have also derived an 
approximation of the spectrum analytically. The results from both calculations are consistent 
with each other. Discussions are given on the possible detections. 

Then I will mention sketchily our analytic calculation of the CMB polarization produced 
by Thompson scattering in the presence of the RGW. The resulting spectra are quite close to 
the numerical one computed from the CMBFAST code, and have several improvements over 
the previous analytic results. Moreover, the formulae bear the explicit dependence on such 
important processes, as the decoupling, the inflation, the dark energy, the baryons. 

2 RGW in the Accelerating Universe 

Consider a spatially flat Universe with the Robertson-Walker metric 



where hij is 3 x 3 symmetric, representing the perturbations, r is the conformal time. The 
scalar factor a(r) is given for the following various stages. The initial stage (inflationary) 



ds 2 = a 2 (r)[dr 2 - (Sy + h^dx^x 4 ], 



(1) 



a(r) = l Q I r 



1+/3 



— OO < T < Ti, 



(2) 



where 1 + (3 < 0, and T\ < 0. The special case of (3 = — 2 is the de Sitter expansion of 
inflation. The reheating stage 

a(r)=a z (r-r p ) 1+ ^, n < r < r s , (3) 

allowing a general reheating epoch. [28, 29] The radiation-dominated stage 

a(r) = a e (r - r e ), r s < r < r 2 . (4) 

The matter-dominated stage 

a(r) = a m (r - r m ) 2 , r 2 < r < r B , (5) 

where te is the time when the dark energy density p A is equal to the matter energy density p m . 
The redshift z E at the time r E is given by 1 + z E = (If 1 ) 1 / 3 - If the current values f2 A ~ 0.7 
and VL, m ~ 0.3 are taken, then 1 + z E ~ 1.33. The accelerating stage (up to the present time 
Ttf)[30]-[32] 

a(r) = l H \r - r a \^, t e <t<t h , (6) 

where 7 is a parameter. For the de Sitter acceleration with = 1 and Vt m = 0, one has 
7 = 1.0. We have numerically solved the Friedman equation 

(^j =H 2 (n A + a rn a' s ) (7) 

with a' = da(r)/dT, and have found that the expression of (6) gives a good fitting with 

7 = 1.05 for tt A = 0.7, 7 = 1.06 for fi A = 0.65, 7 = 1.048 for fi A = 0.75, and 7 = 1.042 for 
Vl A = 0.80. [30]" [32] 

There are ten constants in the above expressions of a(r), except (3 and f3 s , that are imposed 
upon as the model parameters. By the continuity conditions of a(r) and a{r)' at the four given 
joining points t±, t s , t 2 , and t e , one can fix only eight constants. The other two constants 
can be fixed by the overall normalization of a and by the observed Hubble constant as the 
expansion rate. Specifically, we put a(r H ) = l H as the normalization, i.e. 

\TH-T a \ = \, (8) 



and the constant l H is fixed by the following calculation 



To completely fix the joining conditions we need to specify the time instants t\, t 2 , t s , and te- 
From the consideration of physics of the Universe, we take the following specifications[30]~[32]: 

a(r H )/a(r E ) = 1.33, a{r E )/a(r 2 ) = 3454, a{r 2 )/a{T s ) = 10 24 , and a^/a^) = 300. The 
physical wavelength A is related to the comoving wave number k by 

A=^I>. (10) 

The wave number corresponding to the present Hubble radius is k H = 2tt(i(th)/Ih — 2n. 

The gravitational wave field is the tensorial portion of hij, which is transverse-traceless 
dih 1 * = 0, 8^ = 0, and the wave equation is 

fy(v^^Mx,r)) = 0. (11) 

For a fixed wave vector k and a fixed polarization state a = + or x , the wave equation reduces 
to 

fcW" + 2 -h£ Y + k 2 ht ] = 0. (12) 

(X 

Since the equation of h^\r) for each polarization a is the same, we denote h^\r) by h^(r) 
in the following. Once the mode function h k (r) is known, the spectrum h(k,r) of RGW is 
given by 

h(k,r) = ^k\h k (r)\, (13) 



"K 

and the spectral energy density Q g (k) of the GW is defined 

n g (k)^^h 2 (k,r H )(J-J , (14) 

which is dimensionless. 

The initial conditions of RGW are taken to be during the inflationary stage. For a given 
wave number k, the corresponding wave crossed over the horizon at a time r it i.e. when the 
wave length was equal to the Hubble radius: Aj = 2na(Ti)/k to l/Hiji). Now the initial 
condition is taken to be 

/ k \ 2+13 

h(k,r l )=Ai [I -j , (15) 

where the constant A is to be fixed by the CMB anisotropics. The power spectrum for the 
primordial scalar perturbations is P 3 (k) oc \h(k,T H )\ 2 , and its spectral index n s is defined as 
P(k) oc k' 11 "' 1 . Thus one reads off the relation n s = 2(3 + 5. The exact de Sitter expansion of 
(3 = — 2 leads to n s = 1, yielding the so-called the scale-invariant primordial spectrum. 



Any calculation of the spectrum of RGW has to fixed the normalization of the amplitude. 
One can use the CMB anisotropics to constrain the amplitude, receiving the contributions 
from both the scalar perturbations and the RGW. The ratio is defined as 

r = P h /P s , (16) 

the value of which has not been observationally fixed up yet. Here the ratio r is taken as 
a parameter. This will determine the overall factor A in (15). Using the observed CMB 
anisotropies[33, 34] AT/T ~ 0.37 x 1(T 5 at I ~ 2, one has 

h(k H , t h ) = 0.37 x l(T 5 r. (17) 

Then the spectrum h(k,r H ) at the present time t h is fixed. 

Writing the mode function h k (r) = // fc (r)/a(r) in Eq. (12), the equation for /^(t) 
becomes 

fi+(tf-?pjvi k = 0. (18) 
For a scale factor of power-law form a(r) oc r a , the general exact solution is 

Mr) = ci(fcr)2 J a _i (kr) + c 2 (kr)h i_ a (kr) , 

where the constant c\ and C2 are to determined by continuity of the function /^(r) and the time 
derivative (/Xfe(r)/a(r))' at the instances joining two consecutive stages. We have analytically 
solved the equation for the various stages, from the inflationary through the accelerating stage. 
The final expressions are lengthy and we do not write down them here. [30]~[32] However, 
the resulting spectrum will plotted for illustration. Taking the ratio r = 0.37 and 7 = 1.05, 
we have plotted the exact spectrum h(k,r H ) in Fig. 1 for three inflationary models with 
P = -1.8,-1.9, and -2.0, and f3 s = 0.598, -0.552, and -0.689, respectively. We also plot 
the spectrum from the numerical calculation in Fig. 2. And Fig. 3 shows the spectra for 
(3 < —2.0. Fig. 4 compares the root mean square spectrum h(v, T H )/y/v of the model 
7 = 1.05 with the sensitivity of LIGO I SRD [35]~[37] in the frequency range v — 10 ~ 10 4 Hz. 

The spectrum h(u,T H ) depends on the dark energy through the parameter 7. In Fig. 
5 we have plotted the spectra in a narrow range of frequencies. It is seen that the amplitude 
in the model 7 = 1.06 is about ~ 50% greater than that in the model 7 = 1.05. That is, in 
the accelerating Universe with £7a = 0.65 the amplitude of relic GW is ~ 50% higher than 
the one with Qa = 0.7. This difference is probably difficult to detect at present. However, in 
principle, it does provide a new way to tell the dark energy fraction f2 A in the Universe. 
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Figure 1: For a fixed 7 = 1.05, the exact spectrum h(v, r#) is plotted for three inflationary 
models of (3 = —1.8, —1.9, —2.0, respectively. 





1 


1 






1 


1 


1 




1 
























— 








































































































dot line p 


= -1.E 


















































! 


< f i>«, :•' 


7 


sc 


)lid line p = 


-1. 


9 










































* j 


■ l\flA t 

« KM 

k 


*!..'«( 












i 

V;l| 




























dash dot dot line 


3 = -2.C 




! 


I A ll ' 




nl m )' 


."' 




1 










1 


1 








•J-l 1 . '. 


■mm 

1 -!i 'i 


! 



-25 -20 -15 -10 -5 5 10 15 20 



iog 10 v(Hz) 



Figure 2: The numerical spectra h(k,Tn) in the accelerating universe for {3 = —1.8, —1.9 
and —2.0, respectively. 




Figure 3: The exact spectra h(k,Tn) in the accelerating universe for (3 = —2.01,-2.02 
and —2.03, respectively. 




Figure 4: For 7 = 1.05 the root mean square spectrum h(v, Th) / \fv is plotted for the mod- 
els of (3 = —1.8, —1.9, —2.02 to compare with the sensitivity curve from S5 of LIGO.[35] 
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Figure 5: The amplitude of h(u, th) for the model 7 = 1.06 is ~ 50% higher than that of 
model 7 = 1.05. 
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Figure 6: The spectral energy density £l g {v) is plotted for the models of /3 — —1.8, 
P = —1.9, and (3 = —2.0. The model (3 = —1.8 is ruled out by BBN constraint. 



Let us examine the spectral energy density Q g (v) and its constraints. Fig. 6 is the the plots 
of Vt g (v) defined in Eq. (14) for 7 = 1.05. If we use the result LIGO third science run[36, 37] 
of the energy density bound for the flat spectrum with f2 < 8.4 x 10~ 4 in the 69 — 156 Hz 
band, then the model j3 = —1.8 is ruled out, but the models j3 < —1.9 survive. However, this 
LIGO constraint is not as stringent as the constraint by the so-called nucleosynthesis bound. 
[38]"[40] 

J n g (u) d(log v) < 0.56 x 10" 5 . (19) 

Note that this is bound on the total GW energy density integrated over all frequencies. The 
integrand function should also have a bound Vt g {v) < 0.56 x 10~ 5 in the interval of frequencies 
5(\ogv) ~ 1. By this constraint it is also seen from Fig. 6 that only the model (3 = —2.0 are 
still robust. 

We have also obtained the following expressions for the analytic approximate spectrum 

/ u \ 2+l3 

h(k,T H )=Al — \ , k<k E ; (20) 

M^,)^(Ay kE < k < kH ; (21) 

M^)^^__i__ kH <k<k 2 ; (22) 

h{k ' T " ) " A {r H ) ijT^r- k ^ k ^ k - (23) 

h ^ T - )KA (t) wrier" (24) 

where the small parameter e = (1+/3)(1— 7)/7- Approximately , 1 , 3+e ~ , 1 , 3 = VL m /VL^. 
This extra factor reflects the effects of acceleration caused by the dark energy. Some of other 
works on RGW can be found in Refs.[41] — [45], and the effects of neutrino free-streaming have 
been recently computed in Ref.[46] — [49]. 

3 CMB Polarization 

At the beginning, I mentioned that the magnetic polarization of CMB gives another way to 
detect RGW. During the era prior to the decoupling in the early Universe, the Thompson 
scattering of anisotropic radiation by free electrons can give rise to the linear polarization 
only, so we only consider the polarized distribution function of photons / = (I h I r ,U) whose 



components are associated with the Stokes parameters: I — Ii + I r and Q — Ii — I r . The 
evolution of the photon distribution function is given by the Boltzmann equation[50] 



-L _|_ ff — = — — — - q(f - J) 

dr dx i dr dv 



(25) 



where n % is the unit vector in the direction (0,0) of photon propagation, q is the differential 
optical depth and has the meaning of scattering rate. The scattering term q(f — J) describes 
the effect of the Thompson scattering by free electrons, and the term — reflects the effect 
of variation of frequency due to the metric perturbations through the Sachs-Wolfe formula 



1 dv 
v dr 



ldhij ^ ■ 
--— J -nn J . 
2 dr 



(26) 



In the presence of perturbations hij, either scalar or tensorial, the distribution function will be 
perturbed and can be written as 



fo 




+fi 



(27) 



where fi represents the perturbed portion, /o(^) is the usual blackbody distribution. The 
tensorial type perturbations h^, representing the RGW, has two independent, + and x, po- 
larization. 

h ij = ht j + K j = h + e+ j + h x 6* j . 

To simplify the Boltzmann equation (25), for the = h + efj polarization, one writes fi in 
the form [8, 13] 



h = | (I-// 2 ) cos 20 1 







2 



/ (1+^)008 20 

-(1 + /j 2 ) cos 20 
\ 4/isin20 



(28) 



where ( oc I\ + I r = I represents the anisotropics of photon distribution, and f3 oc Ii — I r — Q 
represents the polarization of photons. From the Boltzmann equation, upon taking Fourier 
transformation, retaining only the terms linear in h^, and performing the integration over d\i, 
one arrives at a set of two equations for the + polarization, [13, 18, 51] 

rfln/o 



i k + [ikfjt, + = 



3q 



dhi u 



h + [>h> <i\ h = g / 1 dfi' [(l + /i' 2 ) 2 fa - - (l - /i' 2 ) 2 



(29) 



(30) 



where £ fe = Cfc + Pk, the over dot " • " denotes d/dr. For the blackbody spectrum f(i>o) in the 
Rayleigh-Jeans zone one has ^hMa°l ~ \ Jhe equations are the same for the x polarization. 



In the following we simply omit the sub-index k of wavenumber, the GW polarization notation, 
+ or x, since both h + and h x are similar in computations. In general, it is difficult to give the 
exact solution of (3 and £, but once derived, they can be expanded in terms of the Legendre 
functions 

to*) = £(2/ + m = £(2z + mw, 
i i 

with the Legendre components 

Ur) = I £ dnt(T,n)P«jM), A(r) = \ £ dn0(T,n)Pi(n). (31) 
It can be shown that the spectrum for electric type polarization is given by 



16tt 



(I + 2)(l + l)A- 2 + 6(/-l)(/ + 2)A + 



1(1 - l)A+2 



(2Z — 1)(2Z + 1) (2/ + 3)(2/-l) (2/ + 3)(2/ + l) 



k 2 dk, (32) 



the spectrum for magnetic type polarization is given by 



4tt 



(I + 2)A-i G - l)A+i 



2/ + 1 



2/ + 1 



(33) 



As Eq. (29) shows, one needs the time derivative of h(r) to solve for £ and (3. For both 
polarization, +, x, the wave equation of the relic GW has been given in Eq. (12). The initial 
condition is taken to be 

h(r = 0) = h(k), h(r = 0) = 0, (34) 
with the primordial power spectrum 



tit 



(35) 



where A T is the amplitude, A; = 0.05 Mpc^ 1 is the pivot wavenumber, and n T is the the 
tensor spectrum index. Inflationary models generically predicts nr ~ 0, a nearly scale-invariant 
spectrum. The resulting /i(r rf ) is plotted in Fig. 7. 

One solves the ionization equations during the recombination to give the differential optical 
depth q. Then one obtains the visibility function V(t), 



V(t) = q(r)e 



-k(to,t) 



(36) 



which satisfies J ro V(r)dT = 1 and describes the probability that a given photon last scattered 
at time r, where the optical depth function k(t , r) is related to q(r) by q(r) = — rf/t(r , r)/dr. 
[52]~[54] Fig. 8 shows the profile of V(t) by the the numerical result from the CMBFAST, 




Figure 7: The derivative of RGWs h(rj d ) as a function of k. The solid line is the sudden 
transition approximation, the dash line is that of the WKB approximation, which is nearly 
overlapped with the dot line of the numerical result. 

which is sharply peaked around the last scattering. In calculation it is usually fitted by a 
Gaussian form[16][55] 

V(r) = V(r d )e W (-^f), (37) 

where r d is the the decoupling time, and Ar d is the thickness of decoupling. The WMAP 
data[33] gives Ar d /r = 0.00143. Then, taking V(T d )r = 279 in (37) yields a fitting shown 
in Fig. 8, which has large errors, compared with the numerical one. To improve the fitting of 
V(t), we take the following analytic expressions, consisting of two half-Gaussian functions, [56] 



(T " T d ) 



2 s 



V{r) = V{r d )e W y-^ K ^-y (r < r d ); (38) 

V(r) = nr d )exp(-^^j, (r > r d ); (39) 

with Ar dl /r = 0.00110, Ar rf2 /r = 0.00176, and (Ar dl + Ar d2 )/2 = Ar d . Fig. 8 shows that 
the half-Gaussian model fits the numerical one much better than the Gaussian fitting. This 
difference will subsequently cause a variation in the polarization spectra. 

We now look for an approximate and analytic solution of Eqs. (29) and (30). On smaller 
scales the photon diffusion will cause some damping in the anisotropy and polarization. Taking 




Figure 8: The visibility function V(r) around the decoupling. The half-Gaussian model 
improves the Gaussian model by ~ 11.5%. 

care of this effect to the second order of the small parameter 1/g <C 1, by some analysis one 
arrives at the expression for the component of the polarization 

0i(t o ) = —Ar d i l / drV(T)h(T)ji(k(r-T )) ^e-io^e-^M (40) 
10 Jo Ji x 

where k(t) = k(t ,t), x = k(t')/k(t), and dr' = ^Ar^ as an approximation. The integra- 
tion J dr involving V(r), which has a factor of the form e~~ a( - T ~ Td ' >2 . As a stochastic quantity, 
h{r) contains generally a mixture of oscillating modes, such as e lkr and e~ tkr , and so does 
the spherical Bessel function ji(k(r — t )). Thus h(r)ji(k(r — r )) generally contains terms 
oc e - ibk ^- TQ \ where b G [-2, 2]. Using the formula 

f 00 e- ay2 e ibky dy = e - ^? f°° e'^dy, 

J—oo J — oo 

and V(t) in Eq. (38) and (39), the integration is approximated by 

/"TO . I . /"TO 

/ drV(r)h(T) 3l (k(T - r )) « -D^fa^fa - r )) / drV(r). (41) 
Jo 2 Jo 

where 

£>(£) = I[ e - a ( fcA ^) 2 + e - Q (^) a ], (42) 



and a takes values in the range [0, 2], depending on the phase of h(r)ji(k(T — r )). Here we 
take a as a parameter. For the Gaussian visibility function one would have D(k) = e~ a ^ kATd ^ . 
The remaining integrations in $ is 

f T0 , in \ f°° dx -±kMx =1kM f°° i f°° dx — S-kx 10 i 20 t *n\ 

/ drVir) / — e w K W x e w*W = / due ™ K — e ™ KX = — \n—. (43) 
io V ' Ji x Jo Ji x 17 3 v ; 

This number is the outcome from the second order of the tight-coupling limit, differing the 
first order result y In y in Ref. [55]. Finally one obtains 

1 20 

A(r ) = -In -i l Ar d h(r d )ji(k(T d -T ))D(k), (44) 

which contains explicitly the time derivative h{j d ) of RGW. Substituting this back into Eqs. 
(32) and (33) yields the polarization spectra 

1 / 1 2D\ 2 r 

Cf X = — [—In—) j PUk(r d -r ))\h(r d )\ 2 Ar!D\k)k 2 dk, (45) 

where "X" denotes "G" or "C" the type of the CMB polarization. For the electric type 

„ , , (l + 2)(l + l) . , x 6(/-l)(/ + 2) . , , 1(1-1) / \ fl „, 

PGI{X) = (21-1)(21 + 1) 31 - 2{X) - (21-1)(21 + 3) 31{X) + (2l + 3)(2l + lf + ^ (46) 

and for the magnetic type 

PaW = l^-.M - (47) 

To completely determine C* x above, we need the initial amplitude h(r d ) to be fixed through 
Eq. (35) by the initial spectrum Ph{k) = rP s , associated to the scalar spectrum P s by Eq. 
(16). WMAP observation [57] gives the scalar spectrum 

P s (k ) = 2.95 x 1(T 9 A(ko), (48) 

with k = 0.05 Mpc -1 and A(k ) = 0.8. Taking the scale-invariant spectrum with ut = in 
(35), then the amplitude A T in (35) depends on r. 

The polarization spectrum Cf c of magnetic type, calculated from our analytic formulae 
(45) and from the numerical CMB FAST, [58] are shown in Fig. 9. The approximate analytic 
result is quite close to that of the numerical CMBFAST for the first three peaks that are 
observable. 



The location of the peaks: 




Figure 9: The magnetic polarization spectrum Cf with the ratio r — 1. The solid line 
is the numerical result from the CMBFAST. The upper dot line is the result from the 
half-Gaussian visibility function with a = 1.7, the middle dot line is with a = 2, and the 
lower dot line is the Gaussian fitting with a = 2. The half-Gaussian fitting is better than 
the Gaussian one. 




Figure 10: Cf depends weakly on the dark energy. A smaller Q\ yields a higher ampli- 
tude and shifts the peaks to larger scales. 

In (46) and (47) the spherical Bessel function ji(k(r d —T )) is peaked at / ~ k{jQ—T d ) ~ kr 
for I ^> 1. So the peak location of the power spectra are directly determined by 



C?*<x h(r d ) 2 k 2 D\k)\ k=l/ro . (49) 

The factor D{k) has a larger damping at larger I, so the first peak of the power spectrum has 
the highest amplitude. From our analytic solution one has h(r d ) 2 oc [j2(krd)] 2 , which peaks 
at kr d ~ 3. Thus Cf x peaks around 

I ~ kr ~ 3r /r d . (50) 

A lower dark energy component leads to smaller r /r d . For Q\ = 0.65, 0.73 and 0.80, 
respectively, and with fixed Vlb = 0.044, VL dm = 1 — Vt^ — £lb, a numerical calculation yields 
that To/r d ~ 50.1, 51.3 and 53.6, respectively. A smaller will also shift the peaks h{r d ) 
slightly to larger scales. Together, a smaller will shift the peak of Cf x to larger scales, 
as demonstrated in Fig. 10. This suggests a new way to study the cosmic dark energy. The 
baryon component also influence the peak location. A higher baryon density f2 6 makes the 
peak to large scales, as is demonstrated in Fig. 11. 



The height of amplitude: 




Figure 11: The dependence of Cf on fi& in the ACDM universe with — 0.73, Qdm — 
1 — Q\ — Qb, and r = 1. A larger f2& yields a lower amplitude and shifts the peaks slightly 
to larger scales. 

C xx depend on the decoupling thickness Ar^ and the damping factor D{k): C xx oc 
Ar$D(k) 2 . For a fixed k, the smaller Ar^ leads to a larger D(k). Atj is mainly determined 
by the baryon density f2j of the Universe. A higher corresponds to a smaller At^. The 
total effect is that a higher Q b leads to a lower Cf c , as is shown in Fig. 11. Besides, a smaller 
f2 A yields a higher amplitude, as seen in Fig. 10. 

The influence of the inflation on C xx : 

The inflation models determine the ratio r and the amplitude \h{Td)\ in (35), including the 
spectrum index nx of RGW. A larger r yields a larger \h{jd)\ and a higher polarization. A 
larger n-r yields higher polarization spectra. 

The more recent treatments can be found in Refs. [56], [59] — [62]. 

4 Conclusion and Discussion 

Our calculations of RGW have shown that in the low frequency range the peak of spectrum 
is now located at a frequency v E ~ (%^) 1 ^ 1> h, where v H is the Hubble frequency, and 
there appears a new segment of spectrum between ue and vh- In all other intervals of 
frequencies > v H , the spectral amplitude acquires an extra factor jf 1 , due to the current 



acceleration, otherwise the shape of spectrum is similar to that in the decelerating models. 
The amplitude for the model fi A = 0.65 is ~ 50% greater than that of the model f2 A = 0.7. 
The spectrum sensitively depends on the inflationary models, and a larger j3 yields a flatter 
spectrum, producing more power. Both the LIGO bound and the nucleosynthesis bound point 
out that the inflationary model (3 = — 1.8 is ruled out, but the model (3 = —2.0 is still alive. 
Our analytic polarization spectra of CMB has the following improvements. 

(i) The analytic result of CMB polarization is quite close to the numerical result from the 
CMBFAST code. The dependence of polarization on the dark energy and the baryons are 
analyzed. A smaller f2 A yields a higher amplitude and shifts the peaks to large scales. A larger 
fl b yields a lower amplitude and shifts the peaks to large scales. 

(ii) Our half-Gaussian approximation of the visibility function fits analytically better than 
the simple Gaussian fitting, and its time integration yields a parameter-dependent damping 
factor. This improves the spectrum ~ 30% around the second and third peaks. 

(iii) The second order of tight coupling limit reduces the amplitude of spectra by ~ 58%, 
comparing with the first order. 

(iv) A larger value of the spectrum index n T of RGW and a larger ratio r yield higher 
polarization spectra. 
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